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1 Introduction 



In the article Q new multiple integral representations for the correlation functions of the XX Z 
spin-i Heisenberg chain have been obtained. In the present article, we apply the results of [||] 
to compute the correlation functions of the spin-^ Heisenberg chain at the free fermion point 
A = 0. 

Generically, the Hamiltonian of the finite cyclic XX Z chain with M sites (where M is 
supposed to be even) has the form 

M . h \ 

Hxxz = ^ ( O-^O-^+l + + ^(c^m-^m+l " 1) " l^^^m j • (1-1) 

m=l ^ ^ 

Here, am^'^ denote the local spin operators (Pauli matrices) associated with the m-th site of 
the chain, A is the anisotropy parameter, and h an external classical magnetic field. The 
particularization of this model to the case A = is known as the XX chain (isotropic XY 
model 0): 

A/ / h \ 

Hxx=Y.[ «+i + «+i - 2 ^™ ) • (^-^^ 

m=l ^ ^ 

In spite of the fact that the XX spin-^ chain is equivalent to a model of free fermions, 
its correlation functions are quite non-trivial. They had been studied for a long period by 
numerous authors |^ ^, |6|, 0, ^ and the key results in this field are presently known. It is 
worth mentioning however that the methods used in the works listed above rely essentially on 
the free fermion features of the XX model. Therefore, they cannot be applied to the more 
general XX Z case, at least without significant modifications. On the contrary, we expect our 
present approach, which relies on multiple integral representations of correlation functions, to 
be instructive for the study of the general case as well. 

In 1992 1^], multiple integral representations for the correlation functions of the XX Z chain 
at zero temperature, A > 1 and h = have been obtained from the vertex operator approach. 
Later, in 1996 |l^ (see also ||Tl[] ), similar answers were formulated for the case |A| < 1. A proof 
of these formulas, together with their extension to non-zero magnetic field, has been obtained 



in 1999 13] for both regimes using algebraic Bethe ansatz and the actual resolution of the 



quantum inverse scattering problem [13, It results from these articles that, starting from 
the so-called elementary blocks, one can in principle obtain a multiple integral representation 
for any n-point correlation function of the XX Z chain. 

More precisely, if denotes the ground state, and E'^''^™ the elementary operators acting 
on the quantum space TCm at site m as the 2x2 matrices -Ef^'*^ = the elementary blocks 

for correlation functions are defined as 



m 



2 



The methods developed in |^]-[|T2|] enable one to obtain the quantity ( |1.3| ) in terms of an 
integral with m integrations. It is clear that an arbitrary correlation function in the ground 
state can be expressed as a linear combination of the elementary blocks ( |1.3| ) and, hence, as 
a linear combination of such multiple integrals. It should be stressed however that, although 
these formulas are quite explicit, the actual analytic computation of these multiple integrals 
is missing up to now. Moreover, the evaluation of correlations of physical relevance, like spin- 
spin correlation functions, is a priori quite involved. Indeed, if we consider for example the 
correlation function (afa^), the identity 



m— 1 



(V'.kf^lV'p) = ^El' - Ef) n {Ef + Ef){E}:: - E^2)\^g) (1.4) 

i=2 

shows that the corresponding linear combination of elementary blocks is actually given as a sum 
of 2"* terms. This means that the number of terms to sum up growths exponentially with m, 
which in particular makes it extremely difficult to solve the problem of asymptotic behavior at 
large distance. 

Thus, up to recently, the situation in this field was as follows. On one hand, the free 
fermion limit (A = 0) of the XX Z model was well studied, but the extension of these results to 
the general case came up against serious problems. On the other hand, the multiple integrals 
approach formally provided the possibility to compute the correlation functions for arbitrary 
A; however, because of the technical reasons mentioned above, no result has up to now been 
reproduced via this method, even for the simplest case of free fermiong|^. 

The goal of this paper is to study the correlation functions of the XX chain using the new 
multiple integral representations obtained in [||. In fact, these new representations are nothing 
but re-summations of the multiple integral expressions for the elementary blocks. In particular. 



they enable us to present the spin-spin correlation functions of the type (1^) as a sum of only 
m terms instead of 2"^. We would like to point out that in this paper we do not obtain new 
results, but only reproduce the known answers via a new method. Moreover, we consider here 
the XX model only as a test for the relevance of the formulas obtained in [^. We hope that 
some of the methods developed in the present publication can be applied (perhaps after certain 
modifications) to the general XXZ chain as well. 

This article is organized as follows. In the next section, we introduce some useful notations 
and give the list of formulas obtained in Q for the correlation functions of the XXZ model. In 
Section ^, we compute the correlation function of the third components of spin. The emptiness 
formation probability is considered in Section |^. In Section |5[ we obtain a Fredholm determinant 
representation for the correlation function The asymptotic analysis of this Fredholm 

determinant is performed in Section Some perspectives are discussed in the conclusion. 

^Recently, in ^j, the probability to find in the ground state a string of particles with spin down (the emptiness 
formation probability) was computed at A = by the method of multiple integrals. Let us stress that, contrary 
to the spin-spin correlation functions, this quantity can be expressed as a single elementary block. We consider 
the emptiness formation probability in Section ^ of the present article. 
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2 Correlation functions of the XX Z chain 



For the reader's convenience, we gather in this section the Hst of results obtained in for 
the correlation functions of the XXZ model. Since eventually we study the limit A = 0, we 
hereafter restrict ourselves to the case lAI < 1. 



Let us first of all recall that the Hamiltonian (IJ.) possesses the symmetries 



M/2 

UH^xziA, h)U~^ = -H^^,i-A, h), U=ll (t|^, 



m=l 
M 



(2.1) 



VH^^,{A, h)V-' = H^^,{A, -h), ^ = n 



m=l 



Due to this freedom, there is no common definition of the XXZ model. It means that the 
expressions of correlation functions obtained in different publications may coincide up to a 
common sign and/or the sign of A and h. 

The standard parameterization of the anisotropy parameter is A = coshr/. In the regime 
|A| < 1 the parameter r] is imaginary, and we set rj = —iC, C > 0- The free fermion point A = 
corresponds to ( = Tr/2. 

The general structure of the expressions obtained in |l[| for the spin-spin correlation functions 
is the following: 

K^fm) = E fd-^'' I / ({A,M,4) Sh{{X,z}). (2.2) 

Here a, /3 = x, y, z, and the functions r"^({A, /i, z\) and /({A, z}) are purely algebraic quantities, 
which in particular do not depend on the regime nor on the magnetic field. Their explicit forms 
for specific correlation functions are given below. 

The integration contour surrounds the point Zj = — «C/2 (zj = —itt/A for A = 0), where 
the function /({A,z}) has a pole. All other singularities of the integrand ( |2.2D are outside the 
contour C^. 

The contours C\ and depend on A and h. In all the examples considered below we have 
C\ = C^. For |A| < 1 and h>0, the contour Cx is an interval [—A, A] of the real axis, where 
the value of A is uniquely defined by A and h, although in the general case the dependency 
A = A(A, h) is rather implicit. At A = 0, however, the integration domain can be found 
explicitly from the fact that cosh 2 A = 4/h. Note that if /i — > 0, one has A ^ oo. On the other 
hand, if h approaches its critical value he = 4, then A ^ and all the correlation functions 
become trivial, which comes physically from the fact that the ground state of the Hamiltonian 
(|1.2|) is then purely ferromagnetic. Therefore we consider below only the case < h < he- Due 
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to the symmetry ( p. 11 ), we also do not need to consider negative magnetic field, although all our 
results remain valid for /i < as wellQ. 

Finally, the integrand ( ^ ) contains a function Sh{{\, z}), which also depends on the value 
of the magnetic field. This function is equal to the determinant of a matrix of elements p{Xj, Zk), 
where p{X,z) is the so-called 'inhomogeneous density', solution of the integral equation 

A 

-2Trip{X,z) + j K{X- p)p{p,z)dp = t{X,z), (2.3) 

-A 

with 

zsin2C \ — isinC , , 

smh(A + iQ) smh(A — i(,) sinh(A — zj sinh(A — z — iQ) 

Note that, at z = —i(/2, the function p{X,z) coincides with the spectral density of the ground 
state. In the free fermion limit A = (C = one has K{X) = 0, and thus 

n(A, z) = —t(X, z) = -. (2.5) 

'^^ ' ' 27r ^ ' ^ 7rsinh2(A-z) ^ ' 

After this general setting, let us now be more specific and present the explicit formulas for 
some of the correlation functions of the XX Z chain in the domain |A| < 1. We consider below 
essentially three different cases: 

a) the correlation function = {(yfa^j^i) in a magnetic field; 

b) the emptiness formation probability r(m) in a magnetic field; 

c) the correlation functions g'^ = in zero magnetic field. 

The reader can find the derivation of the corresponding multiple integral representations in . 

a) The correlation function of the third components of spin can be evaluated from the 
generating functional (exp(/5(5i,m))i where Qi^m = \ ~ '^l)> ^ 



{<yK.+i) = ( 2^^m ^ - 4^?m^ + 1 ) (exp(/3Qi,™)) 



q2 q 

/ /3=0 

(2.6) 

Here, the symbols and 2?^ denote respectively the first and the second lattice derivative, 

Vraf{m) = f{m+l)- f{m), vlf{m) = f{m+l)+f{m-l)-2f{m). (2.7) 
The expectation value of the functional (exp(/3Qi^m)) is given by (5.8) of |l|: 

1 / TT ^-^i /" jn X TJ ( sinh(2;a - ^) sinh(Aa + y ) \ 



(exp(/?Q,„))=^-i,/ns/^"A.n 

n=0 ^ '' j=l _^ a=l 



xW„({A}|{z}) •det„ Mjk{{X}\{z}) 



det„. 



p{Xj,Zk) 



(2. 



^Actually the restriction /i > is not necessary. In the case h < 0, the integration contour C\ becomes 
[— A+ ^, -00 + ^] URU [+CX3+ ^, A+ where A is the real positive solution of the equation cosh2A = 4/\h\. 
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Here and further we use the notation det„ for the determinants of nxn matrices. The function 
Wn is defined by 



sinh(AQ - Zb- iQ sinh(zfc - Ag - Q 
^Jl ^J^ sinh(Aa - Afe - iC) sinh(za - Zb - iC) ' 



(2.9) 



and the entries of the matrix Mj^ are 



MjkiW\{z}) = t{zk,Xj) + e^t{Xj,Zk) Yl 



sinh(Aa — Aj — iC) sinh(Aj — Za — iC) 
sinh(Aj — Xa — iC) sinh(za — Xj — i() ' 



(2.10) 



where the functions t{X,z) and p{X,z) are defined in (2.4) and (2.3) respectively. 

For /i = 0, it is more convenient to derive the correlation function (cf ct^_,_]^) from the 
generating functional (exp(/3(5i,m)<7m+i)- 



d 

(0"lf^m+l) = -2^?m-l^(exp(/3Ql,m)f^m+i; 



h = 0. 



(2.11) 



/3=0 



The expectation value of the functional (exp(/?Qi^m)c^+i) is given by (6.7) of Q: 



(exp(/3gi,„)o-f„+i) = 



n=0 ^ ^ 



■j-j- dzj 



27ri 



^n+l^ . ' ^^^^^^^ ~ 2 ) Smh(Aa + 



2 ' 



n 



sinh(Aa + ^) 



-A- sinh(An+i - Zg) _ -A- sinh(A„+i - Za - O 



^^^^ \^sinh(2;a + ^) sinh(Aa - y)/ 

Wn{{X}\U}) 



\sinh(za + f ) / Viri sinh(A„+i - A^) sinh(An+i - A^ - iC) 



xdet^ 



Mjki{X}\{z}) 



det 



n+l 



P{Xj,Zi),. . . , p{Xj ,Zn),p{Xj,-'-^) 



(2.12) 



b) The emptiness formation probability (the probability to find in the ground state a string 
of particles with spin down in the first m sites) is defined as 



1 z 

rim) = ill 

k=l 



(2.13) 



The multiple integral representation of this correlation function can be easily obtained from the 
generating functional (exp(/3(5i,m)) (see equation (C.9) of ||l|): 



r(m)= hm 1- [ -MiM^lL^^t^ p^x„^k) d^'X, 

a<b 



(2.14) 
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where 



^m({A},m)=nn 



smh(Aa - 6) sinh(Aa - Cfe - O det^ t{Xj,Ck) 



a=l b=l 



sinh(Aa - Afe - <) 



n smh(ea - Cfc) 

a>6 



(2.15) 



c) Let us consider finally the correlation functions = {o'f(7^_^_i) and = {aia^_^_i). 
It is easy to see that g^~{h) = g^{—h). For zero magnetic field, these two correlation functions 
coincide. In this case their multiple integral representation is given by (6.13) of [||: 



■j^ ( sinh(Aa + ^ 



Ysinh(Aa 



1 



sinh(A„+i - A„H 



/n+l \ 

n smh(A„+i - Za- ^C) smh(An+2 - Za) 

a=l 

n 

n sinh(A„+i - A„ - <) sinh(A„+2 - K) 

\a=l J 



Here, 



xWn{{X}, {z}) • det„+iMjfc • det„+2 p{\j,zi), . . . , p{Xj, Zn+i), p{Xj , -'4) 



n n+1 

n n sinh(Aa - Zb- iC) sinh(zb - Xa - <) 

a=l b=l 

n n n+1 n+1 ' 

n n sinh(Aa - Afc - iC) n n sinh(za - Zb - <) 

a=l b=l a=l b=l 



and the entries of the (n + 1) x (n+1) matrix M are 



Mjk = t{zk,Xj) -t{Xj,Zk) 



sinh(Aa — Aj — iC) sinh(Aj — Zb — iC) 
^J^ sinh(Aj - Xa- <) sinh(z;, - A^ - <) 



' 2 



(2.16) 



(2.17) 



j < n, 
j = n + 1. 



(2.18) 



To conclude this section, let us recall once more that all the multiple integral representations 
given above hold for arbitrary — 1 < A < 1. In order to particularize these expressions to the 
case of the XX model, one has to set ( = tt/2, A = arccosh(4//i)/2, and to use the expression 
(|2.5D for the inhomogeneous density p{X,z). 



3 Correlation function of the third components of spin 

Let us begin our calculations with the correlation function (o"f(T^_,_x)' which is the simplest 
spin-spin correlation function of the XX model. 
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Observe that, at ( = 7r/2, the equation ( |2.8| ) simphfies drastically. First of all, the matrix 
Mjfc (2.10) becomes proportional to the Cauchy matrix 

'^•MWm) = Jll]'',,y C = |. (3.1) 
and, hence, to compute its determinant one can use the identity 

n 

n smh{xj - Xk) sinh(yfc - Vj) 

<^^^n-^ ^ = — n • (3-2) 

smh(x,- — Vk) -Tj ■ T / \ 

^ ^ ^ n smh{xj - Vk) 

j,k=i 

However, it is more important to notice that det„Mjfc is proportional to (e^ — 1)": this means 
that, if one takes the first (respectively the second) derivative with respect to /? and sets /3 = 
as in (^]^), only the terms n < 1 (respectively n < 2) in the sum ( ^.8| ) do not vanish. Thus, 
after some simple computation, one obtains 

A 

(*,„.) = ji,/.i./^»(.).-"(A)-5^, (3.3) 



A 



{QlJ = iQi,^) + f I d2Af[(^™(z,)^— (A,)) (det: 



-A 

where we have introduced the notation 



sinh(Aj — Zk] 



(3.4) 



sinh(z - f ^ 



smh(2; + ^) 



Let us first consider the integral (3^). As the contour Cz surrounds only the singularity 
z = — z7r/4, where ip{z) admits a pole of order m, the value of the z-integral in ( f3.3| ) is given 
by the corresponding residue at z = — i7r/4. However, this way to compute the 2;-integral is 
not very convenient, especially for large m. Instead, we suggest to deform the original contour 
Cz into an infinite horizontal strip of boundary F given hy Qz = zq — tt and Qz = zq, where 
< zo < 37r/4. Then, obviously. 



ip'^{z)dz r <f"'{z)dz ^ (/9™(z) 

— 2tti Res 



sinh (X — z) J sinh (A — z) sinh (A 

r 



(3.6) 

2 = A 



Since the integrand is a periodic function with period ivr, and since it vanishes at z ^ ±00, it 
is clear that the integral with respect to the new contour F is equal to zero. Thus, to compute 
the z-integral in (|3.3| ), it is enough to take the residue in the second order pole at z = A. The 
remaining integral with respect to A is then trivially computable, and we obtain 

TTl 

(Qim) = — arctan(sinh 2A). (3.7) 

' TT 
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In the next sections, we shall deal with the change of integration variables cosh2A = {cosp)~^. 
Therefore, let us at this stage introduce po such that cosh2A = (cospo)^^- Then (^]^) takes 
the form 

{Ql,m) = (3.8) 
TT 

where po = arccos (^). 

The integral (^|^) can be taken in the same manner. The integration with respect to zi and 
Z2 leads to 



' vr V vr / 4tt^ J \ smh(Ai - A2) / 

-A ^ ^ 

In fact, we do not need to compute {Qf ^) itself, but only its second lattice derivative. Differ- 
entiating ( |3.9D with respect to m, we immediately arrive at 



l^liQlJ = 2 (v)' - ^2^(1 - cos2mpo). (3.10) 
Combining ( |3.^ ) and ( |3.10| ), we finally obtain 

=(^-1) - -^{l - cos2mpo). (3.11) 

It is worth mentioning that, in spite of the fact that we have formally restricted ourselves 
to the case /i > 0, the result ( |3.11| ) remains valid for /i < as well. For zero magnetic field, 
Pq = 7r/2, and the constant contribution to the correlation function disappears. In order to get 
rid of this constant term from the very beginning, it is more convenient to derive (o"f (T^_|_x) from 
the generating functional (exp(/3Qi^m)f^+i) (see ( |2.12| )). Then, for A = 0, the corresponding 
sum reduces to the single term n = 1, which gives 

{^i^m+i) -^f'^'J cosh 2A2 cosh 2Ai ■ sinh2(A2-z) " ^-^'^^^ 
Using the method of calculations described above, we find 

= ^^((-l)'"-l), at h = 0. (3.13) 

4 Emptiness formation probability 



The emptiness formation probability ( p. 13 ) constitute one of the simplest example of correlation 



functions. In particular, unlike the spin-spin correlation functions studied in Sections |3|, ^ and 
^, it can be directly expressed as a single elementary block of the form (|L^), therefore as a 
single (multiple) integral which can be written in the symmetric form (|2.14D . In this section. 
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we explain how to compute this integral in the case A = 0, and how to analyze its asymptotic 
behavior in the limit m ^ oo. 

Setting C = vr/2 in the equations ( pH] ), ( p^ ), we have 

rim)= lim -1 / f-^WJ^}) f J ) ,r.,^ 



a<f) 



and 



Zn^m, m = n n x^'? n - ^^^^^^^^n . (4.2) 

a=i 6=1 cosh(Aa - Afo) ^^^^ V smh 2(Aj - ^k) ) 



Once again, we have to deal with determinants of Cauchy matrices. Computing them via ( |3.2D 
and setting = —iir/A, we obtain 

m 

j n sinh2(Aa - Xb) 
rim) = -jj:—- / d"^A. (4.3) 



m\{2TTy 



-A 



n cosh™(2Aa) 



a=l 



The representation (4.3) can be reduced to a Toeplitz determinant. Indeed, the change of 
variables cosh2Aj = cos^^pj leads to 



PO ^ PQ 

^iri'—m 

Tim) 



m!(27r)" 



/ n ^^V^ ^> = ^;^n I A(e--) A(e^^) (4.4) 



-PO '^>^ -PO 



where A(e^*^) denote Van-der-Monde determinants of variables e^*^-? . Due to the symmetry of 
the integrand with respect to all pj, one can replace one of these Van-der-Monde determinants 
with the product of its diagonal elements multiplied by ml, which gives us 

r{m) = j n e-^('=-i)^"=det^ (^e^(i-i)p.) d^p = det„ ^ j e'^^-^^P dp . (4.5) 

-po \ -po / 



The representation ( [4. 5]) of T(m) as a Toeplitz determinant has already been obtained in 
from the multiple integral representation given in [^]. 

Thus, ( |4.5D provides an explicit expression of the emptiness formation probability, at least 
if m is small enough. However, it is more important to be able to extract the asymptotic 
behavior of r(m) at m — > oo. There exist several ways to do this. Firstly, one can analyze the 
determinant ( |4.5| ) as in . Secondly, the determinant ( |4.5| ) can be transformed to a Predholm 
determinant of a linear integral operator 0, the asymptotic behavior of which can be evaluated 
from the matrix Riemann-Hilbert problem Here we propose a third approach, based on 
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the application of the saddle point method directly to ( |4.3| ). It is possible that this method can 
be used also for the general XX Z model. 



Let us rewrite (4.3) in the following way: 



rim) = ^ I e-^^«^» dX, (4.6) 

D 



where 

m ^ m 

^({A}) = — Vlogsinh2(Aa-Afe)-- Vlogcosh2Aa, (4.7) 



— ' m 

a>b a=l 



and the domain D is defined by —A < Ai < . . . < Xm < A. The integrand in ( [4. 61 ) is positive 
within the domain D and vanishes on its boundary. Moreover, it is not difficult to check that 
the matrix of the second derivatives d"^ S / dXjdXj. is negatively defined. Hence, the integrand 
has a unique maximum in D, which is given by the system 



dS 2 



m— = — V coth(A,- - Aa) - 2 tanh 2A,- = 0. (4.8) 



dXi m 



Following the standard arguments of the saddle point method, we assume that, in the limit 
m oo, the solutions of the system (|4.g| ) are distributed on the interval [—A, A] according to a 
certain density Pq{X). Then, in this limit, ( [4.8D becomes an integral equation for this density: 



tanh 2A = V.P. J coth(A - )u)po(/^) dfi. {4S 



In its turn, the integral (4.6) can be approximated by the value of the integrand in the saddle 
point: 

r(m) ^ 2™%™'^°, m^oo, (4.10) 

where 

A A 

5o = ^ y log sinh^ {X — fi) po{X) po{n) dXdfi — J log cosh 2A • /3o(A) dA. (4.11) 
-A -A 

The analytic expression of the function po{X) can be determined as follows. Setting x = e^^, 
we transform ( |4.9| ) into 

V.P. [' ^My), (4.12) 



a; + 1 Ja x-y 

where po{x) = ^e'^^po{X) and a = e~^^, b = e^^. The solution of the singular integral 
equation ( [4. 12] ) can be obtained in a standard way via the scalar Riemann-Hilbert problem. 
Let us define 
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Then we have 

f6 



27ripo{x) = - / po{x) dx = l. (4.14) 

At the same time, /(x) satisfies the relation 



and, hence, 



2r 

f_{x) + U{x) = ^^, (4.15) 

fix) = fix) |c + 1 r ^ • rYTlVfOT^] ' (^-16) 
I T^iJa y-x iy^ + i)riiy)} 

( \-l/2 

where / (x) = ( (x — a)(x — 6)1 and C is a constant. Substituting this expression into 
(4.14), we eventually obtain 



^°^"^ = ^^V 2(x-a)(b-x) ' 



which results into 



, , 1 cosh A / cosh2A 
^ " ^cosh2AV sinh(A - A) sinh(A + A) ' ^ ' ^ 



This enables us to obtain the analytic expression of in terms of A: 



, / Vcosh 2A \ 

Taking into account that cosh2A = 4//i, we finally obtain the following asymptotic equivalent 
of r(m) in terms of the magnetic field h: 

/4 - /i\ ~ 

r(m)^(— , m^cx). (4.20) 

Thus, this method provides an alternative derivation of the asymptotic behavior of the 
emptiness formation probability |15]. 



5 Correlation function (cr^cr^^^) as Fredholm determinant 



Unlike the correlations of the third components of spin (see Section |3|) , the correlation function 
('^i^'^m+i) remains non-trivial even at the free fermion point A = 0. In this section, we show 
how to compute it from ( 2.16| ) for zero magnetic field. 

Let us first consider the part of the integrand ( 2.16| ) which depends on A„+i and A„+2- 

I - iQ sinh(A„+2 - ^a) ^ 



det„+2 pi\j,Zk) 



/ n+l 

n sinh(A„+i 

a=l 



sinh(A. 



n+l 



An+2) 



n sinh(An+i - Aa - <) sinh(A„+2 - A^) 

\a=l 
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where one should set 2;n+2 = — in the last column of the determinant of densities. We 
can shift the integration contour for \n+2 by — «C) &nd then replace \n+2 by \n+2 — iC in the 
integrand. This changes the sign of the density function, and we obtain 



det„+2 p{><j,Zk) 
sinh(A„+i - Xn+2 + iC) 



/n+l >^ 
n smh(A„+i - Za- tC) smh(An+2 - Za - iC) 

a=l 

n 

n sinh(A„+i - Aa - <) sinh(An+2 - Xa - iC) 

\a=l ) 



We see that for Q = 11/2 the integrand becomes an antisymmetric function of \n+2 and A„+i and, 
hence, the corresponding integral vanishes. It remains then to take into account the contribution 
of the poles which have been crossed during the shift of the A„+2-contour. It is easy to see that 
we have crossed only one singularity, which corresponds to the pole of the function />(A„+2) 
at An 



+2 



-y. The residue in this point gives 



m— 1 



nUn + 1)! / 27ri 
n=o ^ ^ > j=l 



n+l , 



n+l 



-Q j sinh(z^ - f) \ f sinh(Aa + ^ 



a=l \Smh{Za+'-^ 



_^ ysinh(Aa - f 



/n+l \ 

n sinh(Art+i - Za - iC) sinh(za + ' 

a=l 



sinh(A„+i + 



2 ^ 



2 V 



n sinh(A„+i - Xa- «C) sinh(Aa + ^) 



\a=l 



J 



X W„({A}, {z}) ■ det„+iMjfc • det„+i [p{Xj, Zk)] ■ 



(5.1) 



At this stage, we can again use the fact that M and {p{Xj,Zk)) become Cauchy matrices at 
( = it/2. To compute their determinants, it is convenient to use the following modification of 
(U): 



det,, 



sinh2(a;j - yk) 



n cosh(xj - Xk) cosh(yfc - Vj) 
j>k 

n 

2" n cosh{xj - yk) 
j,k=i 



■ det„ 



sinh(xj - yk) 



(5.2) 



Substituting the corresponding expressions of det„+iMjfc and detn+ip{Xj , Zk) into (|5.lD , we 
arrive at 

n+l 



1 """"^ 1 



2i ^ n!(n + 1)! Vtt 

n=0 



1 



a=l 



sinh(A„+i + f ) 



n+l , n+l 



det 



71+1 



sinh(2fc-Ai) 



\ 



sinh(2fe-A„) 
1 

\ sinh{2fe + f ) ) 



(5.3) 
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As we have seen above, the correlation function (cfo"^,^]^) and the emptiness formation 
probabihty T(m) at A = are completely described by only one or two terms at A = 0. The 
peculiarity of the correlation function {crfa^_^_^) is that, even in the limit of free fermions, all 
the terms of the corresponding series survive. Nevertheless, it is possible to express ( |5.3| ) into a 
more compact form. 

The contour integrals with respect to in (|5.3| ) can be easily computed. Due to the 
symmetry of the integrand with respect to all the variables z^, l<A;<n + l, we can make the 
replacement 

1 \ n+l ^ 

a=l 

Then, inserting for each Zk the factors sinh^^(Afc — z^) and ip"^~^{zk) into the fc-th column of 
the remaining determinant, we can integrate separately each of these columns with respect to 
Zk, using the method described in Section |3|: 



det„H 



^ Vsinh(Aj - Zk] 



1 """^ 1 



n=0 



n+1 



-m+1 



a=l 



detyt+i Ujk 
sinh(A„+i + f)' 



(5.4) 



where 



jk 



99™ ^{zk)dzk 



27ri J sinh(2;fc — Aj) sinh(AA; — Zk) 



1 



(/?™ ^{zk)dzk 



1 



27ri / sinh''(2;fc - Aj) 

Lp'^-'^{zk)dzk 



2i(m — 1 



sinh(Aj — Afc) 
cosh(2Aj) ' 

^"^-HXk) 



lixi J sinh(2;fc + ^) sinh(Afc - Zk) sinh(Afc + 



j <n, j ^ k, 



j <n, j = k, 



n + 1. 



(5.5) 

The sum ( |5.4| ) is very similar to the expansion of a Fredholm determinant of a linear integral 
operator. To make it more clear, let us introduce 



V{X,f,)=i 



ip{X)/^{p) 



m — 1 
2 



m — 1 

2 



(5.6) 



vr sinh(A — ^) 

where V{X, A) is defined by continuity from (|5.6| ). Then the equation (|5.4| ) can be written as 

m— 1 ^ „ 

KW) = E;^ /<i"+iA.det„+iCX 
r).=o ■ ;4 



(5.7) 



14 



with 



^(Aj, A„+i) 



U. 



m— 1 

2ismh(A,+i + f ) 



jk 



m — 1 

^ (Afc) 
7rsinh(Afc + f ) 



ITT ^ 
4 ^ 



i, A; < n, 

/c = n + 1, j < n 
j = n + 1, k < n 
j,k = n+l, 



(5.8) 



^ 27rsmh^(A„+i + 

It is shown in Appendix ^ that (5^) results into the derivative of a Fredholm determinant: 



where / denotes the identity operator and 

R{X,fi) 



— det (/ + y(A,/x) + £i?(A,^)), 



(5.9) 



m — 1 

2 



(5.10) 



sinh(A + f )sinh(// + f ) 
The operator I + V + acts on the real axis, and a is an auxiliary parameter. Since i?(A, fi) 



is a one-dimensional projector, the determinant in (5.9) is a linear function of a. Hence, the 
derivative of the determinant does not depend on a. To compare ( |5.9| ) with the result obtained 
in 1^] one can make the standard change of variables cosh2A = cos~^p, cosh2/i = cos~^g. 
Then, the kernel V and the projector R become 



V{p,q) 



sm 



m—1 



(p - q) 



IT sin h{p — q) 



R(p,q) = (_i)™e^(P+9), 



(5.11) 



where the integral operator acts on the interval [— 7r/2, 7r/2]. Finally, replacing p with —p and 
q with —q, we arrive at 



K^„.+i) = (-ir^det(/ 



sm 



m—1 



{p - q) 



vr sin ^{p — q) 27r 



(5.12) 



This formula coincides with the result of up to the factor (—1)™. The existence of this factor 
is due to the fact that we use a different definition for the Hamiltonian, as it was mentioned in 
the beginning of Section |2|. 

6 Long-distance asymptotics of (c^i^cr^+i) 

The leading asymptotic behavior of the correlation function (cxj^cr^^^) was computed in [^, ^. 
Later, in 0], a Fredholm determinant representation of the dynamic temperature correlation 
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function was obtained for an arbitrary value of the external magnetic field. The determinant 
( 5.12| ) appears to be a particular case of this result. To compute its asymptotic behavior at large 



m, one can use the methods of the matrix Riemann-Hilbert problem which were developed in 
1^] to study the dynamic temperature correlations. However, these methods allow to find the 
asymptotics of the determinant only up to a multiplicative constant, whereas the determinant 
( 5.12| ) can be computed explicitly as a finite product of F-functions. In this section, we present 



the corresponding derivation and reproduce the results of the papers 
Observe first that the kernel V{p,q) ( ^.11 ) is degenerated: 



V{p,q) = - ""T/'' ? = --Y^e^^^-'^^^'^-fl (6.1) 
7rsmi(p -q) it 

Thus, the corresponding Fredholm determinant can be reduced to the determinant of a matrix 
of finite size. Indeed, if a kernel K{p, q) has the form Kij), q) = X^fcLi fk{p)9k{Q)-, then 

det(/ + Kij), q)) = det^{5jk + Mjk), (6.2) 

with 



Mjk= / fjip)gk{p)dp. (6.3) 



c 



Here C is the contour where the operator I + K acts. In our case, C = [— 7r/2, 7r/2], and 

fkip) = --e'P^'-'^l k = l,...,m-l, 



vr 



/™(p) = ^e--f, (6.4) 
<7fc(g) = e-*^(^-T), k = l,...,m. 



Thus, for j < m, 



2 j 0, for j — fc even, 

M,k = -6,k - - j (_i)a-fc-i)/2 . ^ _ ^ ^^-^^ 

whereas the elements of the last line of the matrix M are given by 

a J 0, for A; even, 

= \ (-l)(^-i)/2 . 1 , for k odd. ^^-^^ 

Note that the parameter a enters only the last line. Hence, taking the derivative of the de- 
terminant with respect to a, we need to differentiate only the elements of this line, and we 
obtain 

(f^i^^m+i) = — ^det„(ajfc), (6.7) 
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where 



for j < m, and 



djk 



0, 



for j — k even, 
for j — k odd, 

for k even, 
for k odd, 



(6.8) 



(6.9) 



for j = m. 

Our goal is now to compute detm{ajk)- To do this, let us first reorder the columns and the 
lines of the matrix ajk such that it becomes a 2 x 2 block-matrix. One has to move: 

1. the columns with the number 2k to the position k for k = 1, . . . , [y] 

2. the lines with the number 2j — 1 to the position j for j = 1, . . . , [^j ; 

3. the last line with the number m to the position [^j + 1- 

Here, [^] denotes the integer part of m. After these transformations, we arrive at 

/ 



detm(ajfc) = (-1 



im— 1 



det 



0'2j~l,2k 








0.2j-2,2k-l 



, with ao,2A:-l = am,2k-l- (6.10) 



Hereby the sizes of the blocks are [y] x [y] for a2j-i,2k and [^^^] x [^^^] for a2j-2,2fc-i- 
Observe now that 

a2j-l,2k = a2j-2,2k-l = — 7^ -■ (6-11) 



2j -2k -I 



Hence, we obtain 

detm(ajfc) 



-detr 



1 



)-detrm+il 1 ; 
i^J \2j -2k -I 



(6.12) 



It remains to use the analog of (p.2|) for rational functions 



detr 



1 



Xj - Vk 



n i^j - Xk){yk - Vj) 

j>k 

n ■ 

n i^j - Vk) 
j,k=i 



which gives 



[f J if J 



det^(a,fc) = (-1)—! H H 



1 , J — ^ 

1 j=i 



n n 



1 ii ii ^-.i-.i 

2 fc=l J=l J ^ 2 



(6.13) 



(6.14) 
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After the computation of the products with respect to j, we substitute the result into (6.7) and 
eventuahy obtain 



Thus, we have computed the Fredholm determinant ( ^.121 ) as a finite product of F-functions. 
This form enables one to evaluate the large m asymptotic behavior of the correlation function 
((Tj*'cr~_^j^) in a rather simple way (see Appendix ^). The result reads 



f-iy 

{^1 ^m+i) = ^ exp 



/2?7l 

Conclusion 





dt 


{I 





cosh^ t 



i-iy 

8m^ 



+ 0(m-*) . (6.16) 



To conclude this article, we would like to discuss the significance of our results and the per- 
spectives they open concerning the computation of correlation functions in a more general case. 
Actually, the purpose of the work presented here was double. 

Our first goal was to demonstrate that it was really possible in the free fermion limit to 
compute the spin-spin correlation functions using their multiple integral representations. In- 
deed, the ability of this method to provide effective results, even in the simplest case A = 0, 
has for a long time been seriously under question. Thanks to the new formulas obtained in [|l| , 
which correspond in fact to certain re-summations of the elementary blocks [^^[0], we were 
able here to solve this problem. 

The second goal of this paper concerns the possible application of these new integral repre- 
sentations to the general XX Z model. We hope that some of the technical methods presented 
here are not specific to the free fermion point, and that they can also be successfully adapted 
to study a more general case. 

In particular, the method we used here to compute the z- integrals might be quite efficient 
for the evaluation of the long distance asymptotics of the spin-spin correlation functions. Recall 
that, in this paper, we have deformed the original contour Cz into a horizontal strip T of width 
in. At A = 0, the contribution coming from F vanishes due to the periodicity of the integrand. 
In the general case this property is no longer valid since the density function p(A, z) is no longer 
iTT-periodical (except at C = where n is a positive integer). Nevertheless, it seems that one 
can control the order of the contribution coming from F as m — > oo. Indeed, it is possible 



to choose the strip such that \(p{z)\ < 1 (see (3^)) on its boundaries. Then the factor ip'^{z) 



becomes exponentially small uniformly on an arbitrary finite interval of the new integration 
contour for z. Preliminary estimates show that the integrals over F decrease as some negative 
powers of m as m — > oo. Thus, it is very possible that the leading long distance asymptotics of 
the spin-spin correlation functions are given by the residues of the integrand within the strip F. 
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Finally, we would also like to draw the reader's attention on the method used for the evalua- 
tion of the long-distance asymptotics of the emptiness formation probability. In fact, the saddle 
point approach can be applied directly to the representation ( p. 141 ) in the general case as well. 
One can thus expect that the emptiness formation probability decreases as exp{— c(A, h)m?} as 
m oo. The main obstacle to find the coefficient c(A, h) is related to the asymptotic analysis 
of detm t{Xj,(,k) (recall that for A = this determinant is explicitly computable). The problem 
we mention here coincides one to one with the problem of the computation of the partition func- 
tion of the six-vertex model with domain wall boundary conditions [17|, [|^, which was solved 
for the homogeneous case in [^], [^]. For our purpose, it would be desirable to extend these 
results to the inhomogeneous case as well, which is still an open problem. It seems nevertheless 
that the emptiness formation probability admits a Gaussian behavior. 
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A Fredholm determinant 

Let us first recall the definition of a Fredholm determinant: if an operator I + K acts on an 
interval C as 



[{I + KmX) = ct>{\) + J K{X,fi)cl,{f,)dfi, (A.l) 



c 

then its Fredholm determinant is 



det(/ + K) = y" — / dXi---dXn- det„i^(A,-, A^). (A.2) 
^-^ n\ J 



n=0 p 



The series ( ^.7[ ) has almost the same form. To make it exactly the same, we observe that 
nwkUjk < m — 1, and hence deln+iUjk = as soon as n > m — 1. Therefore the sum in 
can be extended up to infinity. We can then use the identity 



det„+iW;,fc = ^Wn+i,„+i - dein{Wjk + aWj^n+iWn+i,k] 



(A.3) 

a=0 
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which holds for any arbitrary matrix W. The sum (^^) can thus be written as 



n=0 



27rsinh^(A„+i + f^) da 



. m — 1 , , m— 1 

/d'^A.det„(y(A,,A.) + £/y(A,,A„,0^ ^ ^'""^^^ ' ^''^^''""^ 



sinh(A„+i + f )sinh(Afc + f ) 



(A.4) 



0=0 



The series (A.4) is now exactly of the form ( [A.2| ), which means that the correlation function 
i'^t'^m+i) ^^^^ represented as the following Fredholm determinant: 



d 



27rsinh^(z^ + f ) da 



. (A.5) 



0=0 



772 — 1 771 — 1 

X det I / + r (A, + f / v(A, 

2-71 J smh(i/ + x) smh(/i + X 

R 

To reduce this determinant to the form (^]^), we use the following lemma: 
Lemma A.l. Let an integral operator 

I + K{X,fi) + a J K{\,v)y{v)x{^i) dv 
c 

acts on an interval C . Hereby the kernel K{X,fi) and the functions x{X), y{X) are such that the 
Fredholm determinant of this operator exists. Then 

J x{u)y{ij)du - -^-^ det ^/ + i^(A, /u) + a J K{X,u)y{u)x{fi) duj 



0=0 



d_ 
da 



det ( / + A'(A, /u) + ay{X)x{fi) 



(A.6) 



Proof. Assume first that det(/ + K) ^ 0, i.e. that there exists the inverse operator 
I — G = (/ + K)^^. Then one can extract the determinant of the operator I + K: 

J x{i^)y{i') du — det I / + A'(A, ^) + a J K{X,v)y{u)x{ii) dv 



\c 



c 



0=0 



det(/ + A') x{iy)y(u) du - det ^ + aj G{X,iy)y{u)x{fi) di^j 



(A.7) 



0=0 
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The operator G(A, v)y{u)x[^) du is a one-dimensional projector, hence, 

det il + aj G{X,u)y{h')x{fj,) du\ =l + aj G{fj,,h')y{i')x{fi) du dfi. (A. 8) 



Substituting this into ( A. 7 ), we obtain 



J x{u)y{u) du — — \ det il + K{X, fi) + a J K{\,v)y{v)x{ii) dv 
\C J \ C / 

det(/ + ir) / (^5{u - fi) - G{fi,uj^y{iy)x{fi) du du- 



a=0 



(A.9) 



c 



In its turn the last equality is equivalent to 



J x{u)y{u) du — — \ det il + K{X, fj.) + a J K{X,u)y{u)x{fi) du 
KC J \ C / 

det(/ + i^:)^ det I I + a ( (I - G)(X, u)y(u)x(fi) du dA . 



a=0 



(A.IO) 



Finally, transforming the product of the two determinants into a single one, we arrive at ( |A.6D . 

If det(/ + K) = 0, we can consider the modified operator / + ^K, where 7 is some complex. 
The Fredholm determinant det(I + jK) is an entire function of 7 and, hence, it has a finite 
number of isolated zeros in any closed domain of the complex plane. Thus, we can choose 7 
such that det(/ + 7K) / 0, repeat all the transformations described above and continue the 
result to the point 7 = 1. Thus, the lemma is proved. □ 

It remains to observe that the structure of the determinant ( |A.5 ) coincides with the one of 



(|AJ| if we set K{X, n) = V{X, /i) and 

x(A) = y(A) 

Thus, we arrive at (|5.9|). 



ip 2 (A) 



27rsinh(A + f ) 



(A.ll) 



B Asymptotic study of the product of F-functions 

We need to compute the asymptotic behavior of the quantity 



n 



Lr(A:-i)r(A; + i; 



(B.l) 
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where — > oo. Expanding (f)^ into a Tailor series, we obtain 



2EE 



(2n)\ V 2 

ln=l ^ ' ^ 



where 



•2n 



logr(z). 



The sum with respect to k in ( |B.2D can be computed using 

N 

E 

k=l 

Substituting this into 



^ V^'^ (k) = iVV'^'^ (A^ + 1) + s (^'("-i) (A^ + 1) - iJ^'-^^ (1) 
we obtain 

^Ar^(2n"l)(Ar + 1) + (2n - 1) fV'(2n-2)^^ ^ ^) _ ^(2n-2)(-^^ 



(B.2) 



(B.3) 



(B.4) 



^ (2n)! V2 

(B.5) 

This series is absolutely convergent, and we can make an asymptotic estimate of each term 
separately. To do this, let us present (pN tbe form 



6^ = ^logiV + 5i + 52, 



(B.6) 



where 



1 °° 9 1 / 1 \ 2" 

S, = -{N^l^'iN + 1) + V.(iV + 1) + C - logiV) - ^-j ^(2n-2)(i)_ 



Here C = — V'(l) is the Euler constant. The last term in ( |B.(]| ) is 

1 /l^ 



52 = 2^. 



iV^(2n-i)(7v + 1) + (2n - l)V'(^'""^)(iV + 1) 



(2n)! V 2 

n=2 ^ ' ^ 

Recall also the asymptotic expansion for the logarithm of F-function 

i 

^ lo? Z-I-— \ci&{1tt\a 



(B.8) 



1 1 ""^ 

logr(z) = zlogz-2--log2:+-log(27r)+^ 



di2k 



k=\ 



2k{2k - l)z2'=- 



j+Oiz'-'-^), \z\ ^ oo, (B.9) 



where i?2fe are Bernoulli numbers. Using ( [B.9D one can easily see that at — > oo the term Si 
has a finite limit, while 5*2 vanishes. 

To compute the limiting value of Si , one can use for example the integral representation for 
the derivatives of ip{z): 



iP^^-\z) 



oo g— 42^2n— 1 

1 - e-* 



dt, n > 1, 3^(z) > 0. 



(B.IO) 
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Then the series in (B.7) becomes the expansion of the exponent, and we obtain 

1 dt 



4./0 t 



-At 1 



cosh^ t 



N^QO. (B.ll) 



This gives us constant the contribution to i;^Ar. 

In order to find the corrections to this quantity, we need first to take into account the higher 
order corrections to tp'^N + 1) and i/;(A^ + 1) in Si and then to take several terms from the series 

In particular, we find 



1 1 dt 



cosh^ t 



This formula can be directly applied for the computation of the asymptotic behavior of ( |6.15D . 
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